Tunable Plasmon Molecules in Overlapping Nanovoids 
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Coupled and shape-tailored metallic nanoparticles are known to exhibit hybridized plasmon reso- 
nances. This Letter discuss the optical properties of a complementary system formed by overlapped 
nanovoid dimers buried in gold and filled with silica. This is an alternative route for plasmon engi- 
neering that benefits from vanishing radiation losses. Our analysis demonstrates the possibility of 
designing artificial plasmon molecules on the basis of void plasmon hybridization, which allows fine 
mode tuning by varying the overlap between voids. The proposed structures could find application 
to both signal processing through buried optical elements and tunable-plasmon biosensing. 
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Plasmons in metallic particles have recently attracted 
considerable attention as promising candidates to real- 
ize photonic devices, enabling manipulation of light over 
sub-micron distances. The frequency of plasmon reso- 
nances can be tuned at will by nanoengineering the par- 
ticle environment and playing with different geometries 
[H, i, S, i, i, i, 0, HI- The suite of available nontriv- 
ial structures includes nanoshells [l], nanorings [1], 
nanovoids nanorods and exotic polyhedral crys- 
tallites . Plasmon hybridization originating in intra- [l| 
and inter-particle [H] interactions adds up an extra han- 
dle to fine tuning resonant mode energies, particularly in 
dimers near the touching limit 0, . 

In this context, the guiding properties of nanoparticle 
linear arrays have been thoroughly investigated with a 
view to designing optical interconnects based upon plas- 
mon hopping between particles 1^, 11 1. However, plas- 
mon propagation in these systems is found to be severely 
damaged by both dissipation in the metal and radiative 
losses. The latter are minimized by relying on small par- 
ticles (< 50 nm), but there is a tradeoff between both 
types of losses. Ultimately, absorption limits propaga- 
tion to sub- micron distances in small particles [11| . This 
problem can be solved by using instead larger voids or di- 
electric inclusions buried in metal. Actually, the radiative 
damping of hybridized void-like plasmons in nanoshells 
decreases extremely rapidly with increasing shell thick- 
ness The basic unit of such void arrays is the void 
dimer, to which we devote this work. 

In the electrostatic limit, the mode frequencies of the 
particle array u/^ and its complementary void array cjj 
are interconnected through the general, exact relation 
(ci -I- eo)/(ei — Co) = [I4I, where e\ and are the per- 
mittivities of the material inside and outside the spheri- 
cal boundaries, respectively, and Aj are eigenvalues that 
depend exclusively on geometry and not on the actual 
dielectric properties of the materials involved. When the 
media under consideration are air and a Drude metal de- 
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one obtains the sum rule (wp^ -I- (i^J)^ 



ujI There- 
fore, the analysis of electrostatic void arrays is fully con- 
tained in the electrostatic particle array. We will however 
consider relatively large voids yielding longer propagation 
distances. Then, retardation effects become dominant for 
sizes above the skin depth (~ 25 nm at visible and NIR 
frequencies), and therefore a separate analysis becomes 
necessary. 

In this work, we describe the optical properties of plas- 
mon molecules formed in overlapping nanovoids buried 
in metal. Similar to particle dimers [5;], large shifts of 
hybridized modes are observed, controlled by varying 
the degree of overlap, while significant field enhancement 
takes place near the junction region. But unlike particle 
dimers, void dimers do not present a singular transition 
near touching conditions. More precisely, the unphysical 
mode observed in the separate particle dimer 7] , charac- 
terized by a net induced charge sitting in each particle, 
becomes now physical for the voids, although the inter- 
void interaction is drastically reduced in non-overlapping 
systems due to the skin-depth effect. 

We rigorously solve Maxwell's equations in the fre- 
quency domain u) for overlapping void dimers filled with 
silica using the boundary element method (BEM) [l3 |. 
in which surface charges and currents are introduced to 
match the fields across dielectric boundaries, thus result- 
ing in a set of surface-integral equations that is trans- 
formed in turn into a linear- algebra problem upon dis- 
cretization of the integrals. The NIR permittivity of gold 
is modeled for simplicity through E g. with parame- 
ters ujp — 7.9 eV and 77 = 0.09 eV [15'|, while we take 
£1 = 2 for silica. 

Fully buried voids are unaccessible to external radi- 
ation, so light scattering cannot be employed to char- 
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with the analytical expression (see arrows) [12 
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FIG. 1: Plasmon modes in overlapping nanovoids. The pho- 
ton local density of states (LDOS) is calculated inside two 
overlapping silica inclusions in gold (solid curves) and pro- 
jected along directions (a) parallel and (b) perpendicular to 
the dimer axis (see insets) . The inclusions are spheres of ra- 
dius a — 240 nm. The overlapping distance is d = 60 nm. The 
LDOS of a single inclusion is shown as dashed curves. The 
vertical arrows mark the dipole plasmon energy of a single 
inclusion obtained from Eq. ([2]) with 1 = 1. 



acterize them. Instead, we investigate the photonic lo- 
cal density of states (LDOS) in the silica region, defined 
as the combined local intensity of all eigenmodes of the 
system under consideration [l6|. With this definition, 
the LDOS p projected along the direction of an atomic- 
decay transition-dipole d is related to the decay rate as 
r = (Att'^lucP /Tl)p. In practical terms, the LDOS can be 
calculated from the field scattered by the void bound- 
aries, projected along the dipole, according to 



1 



Im{d* • E'''=^'} , 



where the first term on the right hand side is the silica 
LDOS. Here, we consider projection directions parallel 
and perpendicular to the dimer axis and give the LDOS 
normalized to its value in silica. 

We show in Fig.[T]thc LDOS at the center of one of the 
silica-filled voids of a dimer formed by spheres of radius 
a = 240 nm with d=60 nm overlap (solid curves). For 
comparison, the isolated-void LDOS is also represented 
(dashed curves), exhibiting a prominent feature corre- 
sponding to the I = 1 electric Mie resonance in agreement 
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where pj — ^yejcua/c, ji and /i^^ are spherical Bessel and 
Hankel functions, and I is the orbital momentum num- 
ber. The small mismatch between BEM calculations and 
the analytical mode frequency originates in non-resonant 
contributions to the dipole-induced field of the former. 

The Z-plasmon mode of the single void is 2Z-I-1 degener- 
ate, and one can expect partial lifting of this degeneracy 
due to void overlap. This is clearly observed in Fig. [T] 
The resulting modes have m = and m = ±1 azimuthal 
symmetry in Figs. HJ^a) and (b), corresponding to tran- 
sition dipoles oriented parallel and perpendicular to the 
dimer axis, respectively. The latter is doubly degener- 
ate, so that the overall number of 6 modes is preserved 
after inter-void interaction is switched on. The dimer 
modes are hybridized and energy split, but the splitting 
is significantly larger for the transverse low-energy fea- 
ture D. The origin of this can be traced back to strong 
interaction within the inter-void opening, which is made 
possible for co-directional transverse orientations of the 
exciting dipoles and has similar nature as gap modes in 
particle dimers In a related context, this gives rise to 
rim modes localized at the opening of voids partly buried 
in flat metal surfaces 
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This intuition is supported by the induced-field distri- 
butions shown in Fig. [2] for features A-E of Fig. [TJ See 
for instance the large increase in field intensity near the 
void junction region associated to peak D. This is ac- 
companied by significant piling up of induced charge, as 
shown in the graph to the right of the contour plot. How- 
ever, peak E has vanishing intensity in the voids opening 
due its quadrupolar symmetry, and therefore it suffers a 
rather moderate shift. 

These results suggest the possibility to design plasmon 
molecules of tailored frequency and symmetry. For ex- 
ample, a quadrupole plasmon mode can be obtained in 
the dimer near the single-void dipole mode frequency (E 
in Figs. [1] and [2]). More sophisticated plasmon modes 
should also be possible by combining several voids, with 
the dimer acting as a building block. 

In contrast to metallic-particle dimers, in which plas- 
mons interact through a dielectric environment Q, our 
void plasmons are coupled only through the region of 
overlap. Therefore, plasmon mixing critically depends on 
the overlap distance d and can be significantly reduced 
in buried geometries, thus simplifying the understanding 
of the resulting plasmon molecule. In particular, modes 
B, C, and E are reliably described in terms of dipole- 
dipole interaction, whereas mode D, which lies further 
away from the single- void dipole resonance, is dominated 
by the rim-like mode, involving a large degree of mixing 
of higher-order modes. 
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FIG. 2; Plasmon mode maps corresponding to peaks A-E in Fig. [T] The contour plots show the square of the induced electric 
field produced by dipoles as indicated by the double arrows. Schemes of the charge distribution are shown to the left of each 
contour plot. The induced surface charge is shown as well, with negative and positive values represented in blue and red, 
respectively. 



The full d dependence of the void dimers is presented in 
Fig. [31 Several characteristic features are clearly resolved 
in the evolution of the plasmon modes: 

(i) Mode splitting is consistently larger for the trans- 
verse orientation of the dipole, in good agreement with 
the conclusions extracted from Figs. [T] and [2l 

(ii) As anticipated above, the interaction between voids 
switches on smoothly between non-touching and overlap- 
ping configurations, in contrast to the singular transition 
observed for particle dimers 0]. 

(iii) The low-energy dipole modes evolve smoothly to- 
wards the single- void dipole when the overlap is complete 
{d — 2a). An initial red shift after touching is followed 
by a blue shift as d increases. In the rim-like mode [see 
D in Figs. [T] and [21 and low-frequency features in Fig. 
[Sjc)], a maximum red shift is produced at an optimum 
distance resulting from the compromise between increas- 
ing mode volume as d becomes smaller (this lowers the 
kinetic energy of the mode) and increasing inter-void in- 
teraction as d increases (this lowers the mode energy due 
to inter-mode attraction). 

(iv) The high-energy quadrupolar modes march to- 
wards the single-void quadrupole as d increases and the 
quadrupole-mode LDOS becomes zero at the center of a 



single- void. 

(v) A weak transition resonance is observed in the 
vicinity of 1.7 eV [Fig. [3{b)], originating in higher-order 
void multipoles that are the counterpart of the left-most 
rim-like mode. 

The plasmon widths determine how well the modes are 
defined. In this sense, the 1.7 eV resonance just noted is 
a well-defined one. In fact, the width of all our reported 
plasmons arise from metal absorption, so that modes as- 
sociated to fields mainly contained in the silica and char- 
acterized by minimum overlap with the metal will have 
longer lifetimes. This is clear when comparing modes D 
and E in Figs. [IJb) and [51 D is broader because it has 
larger intensity in the metal (see region near the void 
junction). However, much of the field intensity lies in 
the silica, and that is why the mode widths are consid- 
erably smaller than 77 in Eq. (|T]). 

Finally, the location of the probing dipole determines 
the actual degree of coupling to the plasmons. This is il- 
lustrated in Fig.m which shows the evolution of the plas- 
mon LDOS along the dimer axis for features B-E of Fig. 
[H Interestingly, only the low-energy hybrid plasmons (B 
and D), which have the original dipole symmetry, con- 
tribute to the LDOS at the dimer center (thick curves). 




FIG. 3: Evolution of void-dimer plasmon energies with overlapping distance d. The modes are probed through the LDOS at 
the center of one of the voids, projected along the dimer axis in (a) and perpendicular to that axis in (b). The void radius 
is a = 240 nm. The vertical arrows mark the energy of the dipole [l — 1) and quadrupole {I = 2) plasmon modes of a single 
inclusion obtained from Eq. ((2}. 
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FIG. 4: LDOS calculated for two overlapping silica inclu- 
sions in gold as a function of position t along the dimer axis. 
The LDOS is projected along parallel (solid curves) and per- 
pendicular (dashed curves) directions relative to that axis. 
Thick and thin curves correspond to low- and high-energy 
resonances, respectively. The void radius is a = 240 nm and 
the overlapping distance is d = 60 nm. 



The high-energy modes have an overall quadrupolar char- 
acter and can only be excited by dipoles displaced with 
respect to the dimer central region. 

In summary, we have investigated the optical prop- 
erties of two overlapping silica inclusions embedded in 
gold and found hybridized plasmon modes, the mixing 
of which can be effectively controlled by changing the 
degree of geometrical overlap. This opens up the possi- 
bility to realize radiative-free plasmon molecules through 
dielectric inclusions buried in metal. Large field enhance- 
ments are observed on resonance near the overlap region, 
which can find application to surface-enhanced Raman 
spectroscopy and other forms of biosensing. The void 
dimer constitutes the basic building block of buried opti- 
cal circuits for future on-chip interconnects and switching 
elements. In contrast to particle arrays, void arrays will 
benefit from minimum cross-talk across the metal and 
vanishing radiative losses, thus allowing a step forward 
in nanoscale integration. 
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